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We present an all-sky search for gravitational waves in the frequency range 1 to 6 kHz during the first calendar 
year of LIGO's fifth science run. This is the first untriggered LIGO burst analysis to be conducted above 3 kHz. 
We discuss the unique properties of interferometric data in this regime. 161.3 days of triple-coincident data were 
analyzed. No gravitational events above threshold were observed and a frequentist upper limit of 5.4 year -1 
on the rate of strong gravitational wave bursts was placed at a 90% confidence level. Implications for specific 
theoretical models of gravitational wave emission are also discussed. 

PACS numbers: 



I. INTRODUCTION 

LIGO (Laser Interferometer Gravitational Wave Observa- 
tory) njj] is composed of three laser interferometers at two sites 
in the United States of America. The interferometers known 
as HI, with 4 km arms, and H2, with 2 km arms, are colocated 
within the same vacuum system at the Hanford site in Wash- 
ington state. An additional 4-kilometer-long interferometer, 
LI, is located in Louisiana's Livingston Parish. The detectors 
have similar orientation, as far as is possible given the cur- 
vature of the Earth's surface and the constraints of the sites 
on which they were built, in order to be sensitive to the same 
gravitational wave polarizations. The relatively large sepa- 
ration between the two sites (approximately 3000 km) helps 
distinguish an actual gravitational wave appearing in both de- 
tectors from local environmental disturbances, which should 
not have a corresponding signal at the other site. GEO 600, 
a 600 m interferometer located near Hannover Germany, also 
operates as part of the LIGO Scientific Collaboration. 

The fifth science run (S5) of the LIGO interferometers was 
conducted between November 2005 and October 2007. LIGO 
achieved its design sensitivity during this run, roughly a factor 
of 2 improvement in sensitivity over the previous S4 run [2]. 
Additionally, S5 was by far the longest science run and had 
the best duty cycle, collecting a full year of livetime of data 
with all 3 LIGO detectors in science mode. This is an order of 
magnitude greater triple-coincident livetime than all previous 
LIGO science runs combined. The analysis discussed in this 
paper uses data from the first calendar year of S5, covering 
data from 4 November 2005 to 14 November 2006. 

Previous all-sky searches for bursts of gravitational waves 
with LSC instruments have been limited to frequencies be- 
low 3 kHz or less, in the range where the detectors are max- 
imally sensitive J2, [H 0, H, la] . The sensitivity above 1 kHz 
is poorer than at lower frequencies because of the storage 
time limit of the interferometer arms, as demonstrated by the 
strain-equivalent noise spectral density curve for HI, H2 and 
LI shown in Figure Q] Shot noise (random statistical fluctu- 




Frequency [Hz] 



FIG. 1: Characteristic LIGO sensitivity curves from June 2006. Shot 
noise dominates the spectrum at high frequencies. 



ations in the number of photons hitting the photodetector) is 
the dominant source of noise above ^200 Hz. 

Despite the higher noise floor, the interferometers are still 
sensitive enough to merit analysis in the few kilohertz regime 
and there are a number of models which lead to gravitational 
wave emission above 2 kHz. As the sensitivity of gravitational 
wave interferometers continues to improve, it is important to 
explore the full range of data produced by them. LIGO sam- 
ples data at 16384 Hz, in principle allowing analysis up to 
8192 Hz, but the data are not calibrated up to the Nyquist fre- 
quency. Thus, this paper describes an all-sky high frequency 
search for gravitational burst signals using HI, H2 and LI data 
in triple coincidence in the frequency range 1-6 kHz. This 
search complements the all-sky burst search in the 64 Hz-2 
kHz range, described in If20ll . 

This paper is organized as follows: Section II describes the 
theoretical motivation for conducting this search. Section III 
describes the analysis procedure. Section IV discusses gen- 
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eral properties of high frequency data and systematic uncer- 
tainties. Section V discusses detection efficiencies based on 
simulated waveforms. Results are presented in section VI, 
followed by discussion and summary in section VII. 
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II. TRANSIENT SOURCES OF FEW-KHZ 
GRAVITATIONAL WAVES 

A number of specific theoretical models predict transient 
gravitational wave emission in the few-kilohertz range. One 
such potential source of emission is gravitational collapse, 
including core-collapse supernova and long-soft gamma-ray 
burst scenarios J7[] which are predicted to emit gravitational 
waves in a range extending above 1 kHz. In a somewhat 
higher frequency regime are neutron star collapse scenarios 
resulting in rotating black holes iSE!]. 

Another potential class of high frequency gravitational 
wave sources is nonaxisymmetric hypermassive neutron stars 
resulting from neutron star-neutron star mergers. If the equa- 
tion of state is sufficiently stiff, a hypermassive neutron star 
is formed as an intermediate step during the merger of two 
neutron stars before a final collapse to a black hole, whereas 
a softer equation of state leads to prompt formation of a black 
hole. Some models predict gravitational wave emission in the 
2-4 kHz range from this intermediate hypermassive neutron 
star, but in many cases higher frequency emission (6-7 kHz) 
from a promptly formed black hole iTToL [Till . Observation of 
few-kilohertz gravitational wave emission from such systems 
would thus provide information about the equation of state of 
the system being studied. 

Other possible sources of few-kilohertz gravitational wave 
emission include neutron star normal modes (in particular the 
f-mode) ITlll as well as neutron stars undergoing torque-free 
precession as a result of accreting matter from a binary com- 
panion |13j| . Low-mass black hole mergers 11411 . soft gamma 
repeaters 111311 or some scenarios for gravitational emission 
from cosmic string cusps [16] are additional possible sources. 
The majority of predicted high frequency gravitational wave 
signals tend to be of a few cycles duration in most scenarios 
since strong signals tend to lead to strong backreactions and 
hence significant damping. 

While there are specific waveform predictions from many 
of these models (some of which are studied in this analysis) 
these models still have substantial uncertainties and are only 
valid for systems with very specific sets of properties (e.g. 
mass and spin). Thus, as has been done previously for lower 
frequencies in each science run, we use search techniques that 
do not make use of specific waveforms. We require only short 
(<Cl s) duration and substantial signal power in the analysis 
band. 



HI. DATA ANALYSIS 

The process of identifying potential gravitational wave can- 
didate events and separating them from noise fluctuations and 
instrumental glitches takes place in several steps. A schematic 



FIG. 2: A schematic of the analysis pipeline. Triggers, which are 
times when the power in one or more interferometer's readout is in 
excess of the baseline noise, are generated using the QPipeline al- 
gorithm El [H, US. Post-processing includes checking for a cor- 
responding trigger at the other site and clustering remaining trig- 
gers into 1 second periods to avoid multiple triggers from the same 
source. Data quality cuts remove triggers from times with known dis- 
turbances which can contaminate the data with spurious transients of 
mundane origin. The remaining triggers are subjected to auxiliary 
channel vetoes and finally a waveform consistency test is performed 
using CorrPower l23ll . 



outline of the analysis pipeline is shown in Fig. [2] Using 
whitened data, triggers with frequency above 1kHz are iden- 
tified separately at the two LIGO sites using the QPipeline 
algorithm lfl8l,IT^,l20ll . then combined with triggers of consis- 
tent time and frequency at the other site in the post-processing 
stage. The data quality cuts and veto stages remove triggers 
correlated with instrumental and environmental disturbances 
that are known to be not of gravitational wave origin. Re- 
maining triggers are then subjected to a final cut based on the 
consistency of the signal shape in the three interferometers 
Q23D . The analysis procedure is described in greater detail in 
the remainder of this section. 

These procedures were developed using time-shifted data 
produced by sliding the time stamps of Livingston triggers 
relative to Hanford triggers with 100 different time-shifts in 
increments of 5 seconds. Applying multiple time-shifts al- 
lows us to produce a set of independent time-shifted triggers 
with an effective livetime much larger than the actual livetime 
of the analysis. Since 5 seconds is much longer than the light 
travel time between the detectors, even after padding for the 
finite time resolution of our search, no genuine gravitational 
wave signals will be coincident with themselves in the time- 
shifted data streams, allowing us to use this set of time-shifted 
triggers as background data. HI and H2 data streams are not 
shifted relative to each other because their common environ- 
ment is likely to produce temporally correlated non-stationary 
noise, meaning that time-shifts between HI and H2 would not 
accurately represent real background. The analysis was tested 
on a single day of data (December 1 1th, 2005), then extended 
to the entire first calendar year of S5. 

GEO 600, a 600 meter interferometer in Germany, was 
also collecting gravitational wave data during this timeframe. 
However, since the smaller GEO 600 interferometer is sub- 
stantially less sensitive than LIGO, including GEO 600 would 
not have caused a substantial increase in overall sensitiv- 
ity. Also, incorporating an additional interferometer not co- 
aligned with the others would have added substantial compli- 
cations to the analysis, especially since the cross-correlation 
test we perform with CorrPower IT23ll is not designed to an- 
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alyze data from detectors that are misaligned. Thus, for this 
analysis, we used GEO 600 data as a follow-up only, to be 
examined in the case that any event candidates were identi- 
fied using LIGO. Virgo lfl7ll . a 3 km interferometer located in 
Cascina, Italy, was not operating during the period described 
in this paper. Joint analysis of LIGO and Virgo data at high 
frequencies will be described in a future publication. 



A. The QPipeline Algorithm 

The QPipeline algorithm is run on calibrated strain data 
l2ltl to identify triggers. Each trigger is identified by a cen- 
tral time, duration, central frequency, bandwidth and normal- 
ized energy. Any trigger surviving to the end of the pipeline 
described in Figure |2] would be considered as a gravitational 
wave candidate event. However, the vast majority of triggers 
generated by QPipeline are of mundane origin. 

Before searching for triggers, QPipeline whitens the data 
using zero-phase linear predictive filtering IT81, 122I1 . In linear 
predictive filtering, a given sample in a data set is assumed to 
be a linear combination of M previous samples. A modified 
zero-phase whitening filter is constructed by zero-padding the 
initial filter, converting to the frequency domain and correct- 
ing^ for dispersion in order to avoid introducing phase errors 

QPipeline is based on the Q transform, wherein the time se- 
ries s(t) is projected onto complex exponentials with bisquare 
windows, defined by central time r, central frequency /□ and 
quality factor Q (approximately the number of cycles present 
in the waveform). This can be represented by the formula 
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Because it uses a set of generic complex exponentials as a 
template bank, QPipeline thus functions much like a matched 
filter search for waveforms which appear as sinusoidal Gaus- 
sians after the data stream is whitened [18] . This bank of 
templates is tiled logarithmically in Q and frequency, but tiles 
at a given frequency are spaced linearly in time. The templates 
are spaced in such a way that we lose no more than 20% of 
the trigger's normalized energy due to mismatches 5t, 5f and 
SQM- 

The significance of a trigger is expressed in terms of its nor- 
malized energy Z, defined by taking the ratio of the squared 
projection magnitude to the mean squared projection magni- 
tude of other templates at the same Q and frequency: 
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where Shi and Sh2 are the power spectral densities of the two 
interferometers and sh2 (/) and Shi (/) are the frequency do- 
main representation of the strain data coming from HI and 
H2. 

The coherent analysis also defines a null stream, which 
is just the normalized difference between the strain data of HI 
and H2. For lower frequency analyses, if the null H- stream 
value is too large the coherent H + stream is vetoed at the cor- 
responding time 12011 . This is because a signal with consis- 
tent magnitude in both detectors should cancel out to zero, 
so a large null stream value indicates an inconsistent signal 
detected by the two interferometers. However, we do not ap- 
ply this null stream consistency veto in the high frequency 
search and simply take the result of the coherent stream as 
the final QPipeline result for the Hanford site, leaving this 
consistency test as part of the follow-up procedure to vet any 
gravitational wave candidates. This is for two reasons: a.) at 
the time the analysis was designed it was feared that substan- 
tially larger systematic uncertainties in calibration at higher 
frequencies mean that the criterion for what constitutes con- 
sistent behavior between the two Hanford detectors would 
have to have been substantially relaxed, and b.) a smoother, 
less glitchy background population makes this consistency 
test only marginally useful (less than a 1% reduction in the 
clustered coincident background trigger rate) above 1 kHz in 
any case. 

For this analysis, we threshold at a normalized energy Z = 
16 for both sites. Along with CorrPower T (defined in section 
III E) this is one of the variables used to tune the false alarm 
rate of the analysis. In the case of Livingston, Z is simply the 
normalized energy coming out of the Q transform, whereas in 
the case of Hanford, this is the normalized energy coming out 
of the coherent stream. 

While lower frequency data are analyzed at 4096 Hz to save 
on computational costs, this search needs the full LIGO rate 
of 16384 Hz in order to analyze higher frequencies. This 
higher sampling rate required computational tradeoffs rela- 
tive to lower frequency analysis. Specifically, data were an- 
alyzed in blocks of 16 seconds rather than 64 seconds due to 
memory constraints. Additionally, the templates applied cov- 
ered signals with Q from 2.8 to 22.6 rather than extending 
to higher Qs in order to reduce the required processing time. 
This choice of Q range is consistent with theoretical predic- 
tions, since the models under study in this frequency range 
generally predict signals of a few cycles. More detailed infor- 
mation on QPipeline can be found in B2011 and 11811 . 



Z=\X\ 2 /(\X\ 2 ). (3.2) 

A gravitational wave signal would appear identical (in units 
of calibrated strain) in the colocated, co-aligned HI and H2 
detectors at the Hanford site. Therefore, a new coherent data 
stream is formed from the noise-weighted sum of the two data 
streams. Mathematically, this can be expressed as: 



B. Post-Processing of Triggers 

After triggers have been identified at both sites, the two lists 
are combined into one coincident trigger list. In order to form 
a coincident trigger, there must be triggers at both sites which 
have time and frequency values consistent with each other. 
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Specifically, the peak times tjj and tx at the Hanford and 
Livingston sites must satisfy the inequality 

\th~t l \ < max((5r ff ,(5T L )/2 + 20ms (3.4) 

where 5th and 8tl are the durations of the two triggers. The 
time of flight for a gravitational wave traveling directly be- 
tween the detectors is approximately 10 ms, so a 20 ms coin- 
cidence window is somewhat padded to allow for misrecon- 
structions in the central time of the waveform. This is a more 
conservative window choice than that of the corresponding 
QPipeline S5 all-sky burst search at lower frequencies lEoll . 
but the difference in the coincidence window has minimal ef- 
fect on the sensitivity of the analysis. 

Similarly, the central frequencies f n and /q j, of the two 
triggers must satisfy the condition: 

\fo.H - /o,l| < max(5/ ,H,<5/o,i)/2. (3.5) 

where Sfo.n and Sfo.L are the bandwidths of the triggers at 
the two sites. This definition is identical to that used in the 
lower frequency analysis. 

Once this coincident list has been obtained, the coincident 
triggers are clustered in periods of 1 second, taking only the 
trigger with the highest normalized energy, in order to elimi- 
nate multiple triggers from the same feature in the data stream. 
The remaining downselected triggers are referred to as clus- 
tered triggers. 

C. Data Quality Cuts 

Data quality cuts are designed to remove periods of data 
during which there is an unusually high rate of false triggers 
due to known causes. An effective data quality cut should re- 
move a large number of spurious background triggers while 
resulting in a relatively small reduction in the livetime of the 
analysis. These cuts are selected from a predetermined set 
of data quality flags, which identify times in which environ- 
mental monitors suggest a disturbance that might influence the 
gravitational-wave readout. The determination of which data 
quality flags to apply is made based on single detector prop- 
erties and an exact procedure for application of these flags 
is put in place before generating coincidences which may be 
considered gravitational wave candidates. The application of 
data quality flags therefore does not affect the statistical valid- 
ity or "blindness" of the search. We use the same category 1 
and category 2 data quality cuts as the S5 low frequency burst 
searches 112011 . Category 1 cuts remove periods of time where 
there were major, obvious problems, such as a calibration line 
dropout or the presence of hardware injections, which make 
the data unusable. Similarly, category 2 cuts remove periods 
for which there is a clear external disturbance which distorts 
the data. Category 2 cuts result in a loss of 1 .4% of the triple- 
coincident livetime. While category 1 periods are removed 
before the start of the analysis, category 2 periods are removed 
at a later stage so as to avoid creating a large number of very 



short science segments which are impractical to process using 
QPipeline. 

Category 3 data quality flags, which define periods where 
the data are analyzable but still somewhat suspect due to some 
known cause, were studied one at a time for their effectiveness 
relative to high frequency triggers. The category 3 flags used 
for this high frequency analysis are a subset of those adopted 
at low frequencies. Flags which removed QPipeline back- 
ground triggers at a much higher rate than expected by ran- 
dom Poisson coincidence were selected for use. Specifically, 
the rate of clustered single site triggers must be at least 1.7 
times higher for periods when a given data quality flag is on 
relative to periods when that flag is off. As in the lower fre- 
quency analyses, category 3 data quality flags are used only 
for purposes of setting the upper limit, but triggers surviving 
to the end of the pipeline may still be examined as gravita- 
tional wave event candidates if they are within a category 3 
data quality segment. The flags used in this high frequency 
analysis are summarized in Table U Applying the selected 
category 3 data quality flags ultimately removes 19.4% of the 
surviving coincident time-shifted background triggers and re- 
sults in a 1 .7% reduction in triple coincident livetime. 



D. Auxiliary Channel Vetoes 

The LIGO interferometers use a large set of auxiliary de- 
tectors to determine when potential event candidates are the 
result of environmental causes (such as seismic activity or 
electromagnetic interference) or problems with the interfer- 
ometer itself rather than actual gravitational waves. Triggers 
from these auxiliary detectors act as vetoes, removing poten- 
tial gravitational wave candidate events that occur at the same 
time as the trigger in the auxiliary detector. These vetoes are 
distinguished from the data quality cuts described in the pre- 
vious section because they are determined in a statistical way 
and remove triggers from a much shorter period of time (tens 
to hundreds of milliseconds around a particular veto trigger 
rather than blocks of seconds to thousands of seconds in the 
case of data quality cuts). As with the data quality flags de- 
scribed above, all tuning of event-by-event vetoes is done on 
a single instrument basis before coincident triggers are gener- 
ated. Vetoes are divided into categories using the same defini- 
tions as data quality flags. The same list of category 2 vetoes 
used at low frequencies ll20ll was applied to this search. These 
vetoes require multiple magnetometer or seismic channels at 
a given site to be firing simultaneously. 

This analysis also uses the same method of selecting which 
category 3 auxiliary channel vetoes to apply as was used for 
the lower frequency S5 all-sky searches, but used an inde- 
pendent set of high frequency QPipeline time-shifted back- 
ground triggers to select these vetoes. A list of potential 
vetoes is assembled from the various auxiliary channels at 
different thresholds and with different coincidence windows. 
The effectiveness of each potential veto is measured by its 
efficiency-to-deadtime ratio, which is the percentage of back- 
ground triggers it removes from the analysis divided by the 
percentage of the total livetime it removes. The vetoes which 
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TABLE I: Category 3 data quality cuts for high frequency analysis. 



Flag name 


Description 


Livetime 
Loss (s) 


Ratio of Clustered Trigger Rate 
(Flag OmFlag Off) 


H 1 : WIND.OVER.30MPH 


Heavy wind at 
ends of HI arms 


5531 


1.93 


H1:DARM_09.1 1 _DHZ_HIGHTHRESH 


Up-conversion of seismic 
noise at 0.9 to 1.1 Hz 


6574 


1.76 


Hl:SIDECOIL_ETMX_RMS_6HZ 


Saturation of side coil 
current in HI X end mirror 


1360 


2.11 


H 1 :LIGHTDIP_02_PERCENT 


Significant dip in stored 
laser light power in HI 


34336 


2.24 


H2:LIGHTDIP.04_PERCENT 


Significant dip in stored 
laser light power in H2 


40562 


2.04 


L 1 :L1(jH 1 Dlr_04_FbKCbN 1 


Significant dip in stored 
laser light power in LI 


1 15584 


2.85 


L1:BADRANGE_GLITCHINESS 


Abrupt drop in interferometer 
sensitivity, quantified in terms of 
effective range for inspiral signals 


3185 


1.95 


L 1 :HURRIC ANE.GLITCHINESS 


Hurricane was active near Livingston 


42917 


2.92 



are actually applied are selected in a hierarchical fashion, first 
picking the most effective veto, then calculating the effective- 
ness of the remaining possible vetoes after this one has been 
applied. The next most effective veto is then selected and the 
process repeated until all remaining veto candidates have ei- 
ther an efficiency-to-deadtime ratio less than 3 or a probabil- 
ity of their effect resulting from random Poisson coincidence 
greater than 10~ 5 . The vetoes were selected using a set of 
background triggers obtained from 100 time-shifts of LI with 
respect to H1H2, with offsets ranging from -186 to 186 sec- 
onds in increments of 3 seconds. Time-shifts which were also 
divisible by 5 and thus present in the set used to determine 
the final background of the analysis were omitted, making 
the veto training and test sets independent. Of 18831 triggers 
remaining in time-shifted background after category 3 data 
quality cuts, 2284 are removed by vetoes (12% efficiency), 
while the vetoes cause a 2% reduction in the overall livetime 
of the analysis. 



E. Cross-Correlation Test with CorrPower 

The remaining clustered triggers are next subjected to 
cross-correlation consistency tests using the program Corr- 
Power 023I1 . CorrPower has previously been used in S3 and 
S4 analyses HI @]. Unlike QPipeline, which only looks for 
excess power on a site-by-site basis, CorrPower thresholds 
on normalized correlation between data streams in different 
detectors. CorrPower was selected for use in this analysis 
because it is relatively fast computationally and effective for 
roughly co-aligned interferometers such as LIGO. For analy- 
ses including detectors with substantially different alignments 
relative to LIGO, such as Virgo or GEO 600, one does not 



necessarily obtain consistent correlated signals between inter- 
ferometers and more sophisticated fully coherent techniques 
such as Coherent WaveBurst If24ll or X-pipeline If25ll would be 
preferable. 

Before applying the correlation test, data was filtered to the 
1-6 kHz target frequency range of the search. Additionally, 
triggers were rejected entirely if their central frequency as de- 
termined by QPipeline was greater than 6 kHz. Since this 
analysis extends CorrPower to higher frequency regimes com- 
pared to previous analyses, it was necessary to add Q = 400 
notch niters at frequencies of 3727.0, 3733.7, 5470.0 and 
5479.2 Hz, which correspond to "butterfly" and "drumhead" 
resonant frequencies of the interferometers' optical compo- 
nents. The data are whitened. CorrPower then measures cor- 
relation using Pearson's linear correlation statistic: 



x){Vi - y) 



(3.6) 



where x and y are in this case the time series being compared 
for the two interferometers, x and y are the average values and 
N is the number of samples within the window used for the 
calculation. This r-statistic is calculated over windows of du- 
ration 10, 25 and 50 ms. This variable is maximized over var- 
ious time-shifts between the two interferometers. The max- 
imum time-shift between one of the Hanford detectors with 
the detector at Livingston is 11 ms, whereas the maximum 
time-shift between the two Hanford detectors is 1 ms. The 
final output of CorrPower which we use as a data selection 
criterion is called T. T is an average of the r-statistic values 
for each of the 3 detector combinations, using the integration 
length and relative time-shift between interferometers which 
results in the highest overall r-statistic value. 



8 



F. Tuning for the Final Cut 

CorrPower was run on the triggers resulting from the 100 
background time-shifts. This distribution was used to deter- 
mine the value of the cut on the CorrPower T output vari- 
able. In order to obtain an estimated false alarm rate (FAR) of 
around one tenth of an event candidate in the analysis of time- 
shift-free foreground data, cuts were applied to remove the 
bulk of the time-shifted background distribution, only keep- 
ing triggers with V values greater than 6.2 and a Qpipeline 
normalized energy greater than Z = 16 at both sites. This 
results in a final false alarm rate of ~ 10~ 8 Hz. 

IV. PROPERTIES OF LIGO DATA ABOVE 1 KHZ 

A. High Frequency Trigger Distributions 

Although the sensitivity of the detector is poorer at higher 
frequencies, the noise is more stationary in the shot-noise 
dominated regime. QPipeline normalized energy distributions 
from H1H2 for both high (>1 kHz) and low «1 kHz) fre- 
quency triggers are shown for a single day (December 11th, 
2005) in Figure [3] The distribution of single interferome- 
ter triggers at higher frequencies falls off substantially more 
sharply than does the lower frequency distribution and con- 
tains far fewer statistical outliers. The poorer statistics of the 
low frequency data set are due to glitches in the band below 
200 Hz. 



B. Systematic Uncertainties 

Due to variations in the response of the detectors as a func- 
tion of frequency, systematic uncertainties are calculated sep- 
arately for each of three detection bands: below 2 kHz, 2 to 4 
kHz and 4 to 6 kHz. The dominant source of systematic uncer- 
tainties is from the amplitude measurements in the frequency 
domain calibration. The individual amplitude uncertainties 
from each interferometer - of order 10% - are combined into a 
single uncertainty by calculating a combined root-sum-square 
amplitude SNR and propagating the individual uncertainties 
in this equation assuming each error is independent. In addi- 
tion to this primary uncertainty, there is a small uncertainty 
(3.4% or less depending on frequency band) introduced by 
converting from the frequency domain to the time domain 
strain series on which the analysis was actually run 12111 . 

There is also phase uncertainty on the order of a few de- 
grees in each interferometer and in each frequency band, aris- 
ing both from the initial frequency domain calibration and the 
conversion to the time domain. However, phase uncertainties 
are within acceptable tolerance. In this analysis in particular, 
the omission of the null stream in QPipeline means the analy- 
sis is generally insensitive to phase shifts between the interfer- 
ometers on the order of those observed. Likewise, CorrPower 
is mostly insensitive to phase shifts between interferometers 
because it automatically maximizes over multiple time-shifts 
between the interferometers and will therefore still find the 




log 10 (Z) 

FIG. 3: Normalized energy Z of high and low frequency QPipeline 
triggers. The low frequency distribution contains a substantially 
higher number of outliers. 



maximum possible correlation. Some distortion in the shape 
of broadband signals due to differing phase response at differ- 
ent frequencies is in principle possible. However, this is not 
a significant concern since the phase uncertainties at all fre- 
quencies correspond to phase shifts on the order of less half a 
sample duration. We therefore do not make any adjustment to 
the overall systematic uncertainties due to phase error. 

The antenna pattern for LIGO is normally calculated us- 
ing the long wavelength approximation, which assumes the 
period of oscillation of a gravitational wave is large with re- 
spect to the transit time of a photon down the length of the 
interferometer arm and back. This assumption is less accu- 
rate as the frequency increases. However, comparing results 
using the approximate long wavelength antenna pattern and 
frequency-dependent exact antenna pattern lf27tl even towards 
the extreme high end of our frequency range (at 6 kHz) results 
in sensitivity calculations (see next section) differing by only 
^1%. Thus, the approximation of a constant antenna pattern 
has a negligible effect on the analysis. Finally, we include 
a statistical uncertainty of around 2.7% (with some variation 
from waveform to waveform due to different numbers of in- 
jected waveforms). 

In each frequency band the frequency domain amplitude 
uncertainties are added in quadrature with the other smaller 
uncertainties to obtain the total uncertainty. The total 1 a un- 
certainties are then scaled by a factor of 1.28 to obtain the 
factor by which our h ISS limits are rescaled in order to ob- 
tain values consistent with 90% confidence level upper limits. 
These net uncertainty values are 1 1.1% in the less than 2 kHz 
band, 12.8% in the 2-4 kHz band and 17.2% in the 4-6 kHz 
band. Waveforms with significant signal content in multiple 
bands are considered to be in the band with the larger uncer- 
tainty. 
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FIG. 4: Two example high frequency waveforms resulting from grav- 
itational collapse of rotating neutron star models 1 9]. Dl results from 
a nearly spherical 1.26 solar mass star while D4 results from the col- 
lapse of a maximally deformed 1.86 solar mass star into a black hole. 
The figures show the plus polarization for each waveform (the cross 
polarization is at least an order of magnitude weaker in both cases) 
at a distance of 1 kpc, assuming optimal sky location and orienta- 
tion. At this distance, the h ISS magnitudes of the two waveforms are 
5.7 x 10 _22 Hz- 1/2 for Dl and 2.5 x 10" 21 Hz- 1/2 for D4. They 
differ from the figures presented in |9] in that the non-physical con- 
tent at the beginning of the simulations has been removed. 



V. DETECTION EFFICIENCY 

Efficiency curves have been produced for three types of sig- 
nal. The cuts were developed on a set of 15 linearly polarized 
Gaussian-enveloped sine waves (sine-Gaussians) of the form 



h (t +t)=h sin (27r/ t) cxp (- (27r/ i) 2 /2Q' 



(5.1) 



where fo and to are the central frequency and time of the 
waveform and Q is the quality factor defined previously. Ad- 
ditionally, we tested a set of three linearly polarized Gaussian 
waveforms as well as two waveforms taken from simulations 
by Baiotti et al. (9J], which models gravitational wave emis- 
sion from neutron star gravitational collapse and the ringdown 



of the subsequently formed black hole using polytropes de- 
formed by rotation. The two scenarios studied here are des- 
ignated Dl, a nearly spherical 1.26 solar mass star, and D4, a 
1.86 solar mass star that is maximally deformed at the time of 
its collapse into a black hole. These two waveforms are shown 
in Figure H These two specific waveforms represent the ex- 
tremes of the parameter space in mass and spin considered by 
Baiotti et al. 

The BurstMDC and GravEn packages 12811 were used to 
create simulated gravitational-wave "injections" which were 
superimposed on real data in a semi-random way at intervals 
of approximately 100 seconds. This placed all injections far 
enough apart that whitening and noise estimation using data 
surrounding one injection is never affected by a neighboring 
injection. Each waveform was simulated between 1000 and 
1200 times for each of the 18 different amplitudes. The in- 
trinsic amplitude of a gravitational wave at the Earth, without 
folding in antenna response factors, is defined in terms of its 
root-sum-squared strain amplitude: 



h ISS = ^ J (|M*)l 2 + IM0l 2 )dt (5.2) 

where h + (t) and h x (t) are the plus and cross-polarization 
strain functions of the wave. Since h is a dimensionless quan- 
tity, h rss is given in units of Hz -1 / 2 . 

The injections were distributed isotropically over the sky. 
Thus, even a few nominally very strong software injections 
are missed by the pipeline because they are oriented in a very 
sub-optimal way relative to at least one interferometer. Since 
they are simulating an actual astrophysical system, the Dl and 
D4 waveforms also include a randomized source inclination in 
addition to random sky location and polarization. A sin 2 (t) 
dependence on the inclination angle was assumed. Figure 
|5] shows efficiency curves for some of these waveforms as a 
function of signal amplitude. The h ISS values for which 50% 
and 90% of sine-Gaussian injections are detected are summa- 
rized in Table IITI Figure [6] shows the detection efficiency for 
the simulated Dl and D4 Baiotti et al. models as a function 
of distance from Earth, indicating that a neutron star collapse 
would have to happen nearby (within a kiloparsec) to be de- 
tectable at our current sensitivity. 

Hardware injections, wherein actuators were used to phys- 
ically simulate a gravitational wave in the interferometers by 
moving the optical components, were performed throughout 
S5. Although the numbers and variety of amplitudes were 
not sufficient to produce hardware injection efficiency curves, 
sine-Gaussian hardware injections at 1304, 2000 and 3067 
Hz were reliably recovered using the high frequency search 
pipeline at amplitudes large enough that their detection is ex- 
pected based on sensitivities determined by software injection 
efficiencies. Table ITTI1 shows the central frequency, amplitude 
and fraction of hardware injections detected. For hardware in- 
jections, amplitude is given in terms of ft. rss ,det, the root sum 
square of the strain in the detector. This is defined analogously 
to equation |5.2| with /i+.dct and /i x ,det in place of h + and h x . 

Good timing and frequency reconstruction help improve 
detection efficiency. Using Q — 9 sine-Gaussian waveforms, 
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TABLE II: h^T° and h^T" values (the root sum square strain at 
which 50% or 90% of injections are detected) for Q = 9 sine- 
Gaussians. Values in this table are adjusted for systematic uncer- 
tainties as described in section IV B. 



Central Frequency 


h™ % (Hz" 1 / 2 ) 

rss V / 


(Hz" 1 / 2 ) 

rss V / 


1053 


2.87x 10~ 21 


1.97x 10~ 20 


1172 


3.15x 10~ 21 


2.04x 10~ 20 


1304 


3.31 xlO -21 


2.06x 10 -20 


1451 


3.73x 10~ 21 


2.33 x 10~ 20 


1615 


3.99xl0 -21 


2.67 xl0~ 20 


1797 


4.91 xlO" 21 


3.10xl0- 20 


2000 


5.22xl0~ 21 


3.30xl0~ 20 


2226 


6.08xl0~ 21 


3.74xl0~ 20 


2477 


6.63xl0 -21 


4.47 xlO -20 


2756 


7.59xl0~ 21 


5.14xl0~ 20 


3067 


9.20xl0 -21 


5.62xl0~ 20 


3799 


1.17xl0 _20 


8.06xl0~ 20 


3900 


1.19xl0~ 20 


7.87xl0~ 20 


5000 


1.67 xlO" 20 


9.47 xlO -20 



TABLE III: S5 Q = 9 sine-Gaussian hardware injections above 1 
kHz. Note that /i rss ,dct and h ISS are different quantities because h ISS 
does not include a sky location dependent antenna response factor, 
which will reduce the detector response by an additional factor of 
0.38 on average. Care should therefore be taken when comparing to 
Table I. 



Central Frequency (Hz) 


'irss.dct (HZ 1 / 2 ) 


fraction recovered 


1304 


5.00xl0~ 22 


0/2 


1304 


1.28xl0~ 20 


16/16 


1304 


2.56xl0~ 2() 


16/16 


2000 


6.00xl0~ 22 


0/102 


2000 


l.OOxlO -21 


0/4 


2000 


1.20xl0~ 21 


14/127 


2000 


2.40xl0~ 21 


125/125 


2000 


4.80xl0~ 21 


117/117 


2000 


9.60xl0~ 21 


21/21 


2000 


1.92xl0~ 20 


16/16 


2000 


3.84xl0~ 2() 


16/16 


3067 


7.21 xl0~ 21 


13/13 


3067 


1.44xl0~ 20 


13/13 


3067 


2.88xl0~ 2() 


3/3 


3067 


5.76xl0~ 20 


3/3 
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FIG. 5: Software injection efficiency curves for the set of sine- 
Gaussians of various frequencies (top) and Gaussians plus astrophys- 
ical waveforms (bottom). There is a consistent reduction in efficiency 
as a function of frequency following the noise distribution. 
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the timing resolution has been demonstrated to be within one 
cycle of the waveform and frequency resolution is better than 
10%, limited by the coarseness in frequency space of the tem- 
plates used in QPipeline. 



FIG. 6: Efficiency as a function of distance from Earth for supernova 
collapse waveforms Dl and D4 |9], assuming random sky location, 
polarization and inclination angle t. A sin 2 (t) dependence on the 
inclination angle was assumed. 
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FIG. 7: Upper limit curves for a number of our tested waveforms. 
The rate at Earth of gravitational waves of each given type is ex- 
cluded at a 90% confidence level. The curves have been adjusted to 
account for systematic uncertainties as described in section IV B. 



VI. RESULTS 

Having tuned the analysis on background from 100 time- 
shifts and tested it on a single day of data, we then performed 
the analysis on the actual coincident (or "foreground") data. 
No event candidates above our threshold were observed. 

As in previous burst analyses (e.g. J2l), we set single-sided 
frequentist upper limits on the rate of gravitational wave emis- 
sion. The upper limits in the frequency range 1-6 kHz are 
shown in Figure |7]for a sub-sample of our tested waveforms. 
161.3 days of triple-coincident livetime were analyzed (see 
ll29tl for a complete list of analyzed times). After performing 
predetermined category 3 data quality cuts and vetoes, 155.5 
days of triple-coincident data were used to set upper limits on 
gravitational wave emission. 

For gravitational waves with amplitudes such that detection 
efficiency approaches 100%, the upper limit asymptotically 
approaches a value of 0.015 events per day (5.4 events per 
year), as determined primarily by the livetime of the analy- 
sis. While other untriggered searches for gravitational waves 
with comparable or greater livetime (e.g. the corresponding 
LIGO lower frequency analysis Ir20ll and searches by IGEC 
ll3"oll ) have been conducted in overlapping frequency bands, 
this analysis represents the first limit placed on gravitational 
wave emission over much of the frequency band. 

The number of triggers surviving through each stage of the 
analysis are shown in Table [IV] While there are no event 
candidates above our threshold in this analysis, the rates be- 
fore the final CorrPower cut are slightly higher than expected. 
However, assuming Poissonian statistics, this is not a statis- 
tically significant excess since there is a 6.2% chance of get- 
ting at least the observed 193 foreground triggers after all data 
quality cuts and vetoes have been applied. Figure [8] demon- 
strates that the rate of triggers per time-shift can in fact be 
treated as a Poisson distribution. 

The foreground to background consistency of the Corr- 



TABLE IV: Number of triggers surviving various stages of the anal- 
ysis: initial coincident triggers, triggers remaining after the removal 
of segments removed due to data quality criteria, triggers remain- 
ing after vetoes based on auxiliary channels have been applied and 
triggers ultimately surviving after the CorrPower linear correlation 
cut (r). Shown are results for 100 time-shifts, the same result nor- 
malized to the actual livetime, and the foreground results from the 
analysis performed without time-shifting the data. The background 
normalization reflects the fact that the livetime is different for differ- 
ent time-shifts. 





Background 
Count 


Normalized 
Background 


Unshifted 
Count 


Coincident triggers 


23361 


242.9 


265 


After data quality cuts 


18831 


195.8 


223 


After auxiliary channel vetoes 


16547 


172.0 


193 


After T >6.2 threshold 
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FIG. 8: Histogram showing number of time-shifts vs. counts normal- 
ized to analysis livetime. Superimposed is the expected distribution 
based on Poissonian statistics, which is consistent with the observed 
distribution. The black line at 193 counts indicates the actual number 
of foreground triggers observed. 



Power T distribution (Fig. |9]l and QPipeline normalized en- 
ergies from the Livingston and Hanford sites (Fig. ITOb were 
also studied. These plots are produced after all data quality 
cuts and vetoes were applied, but before the final CorrPower 
r cut. The distributions are plotted cumulatively, i.e. each bin 
shows foreground and time-shifted background counts greater 
than or equal to the marked value. Other than the upward 
fluctuation in total counts already discussed, the distributions 
themselves are essentially consistent with expectation. 

Since they appeared to stand out slightly from the expected 
background distribution (although not at a statistically sig- 
nificant level), the loudest 3 triggers in Hanford QPipeline 
normalized energy, loudest 2 triggers in Livingston normal- 
ized energy and the trigger with highest CorrPower T value 
were studied on an individual basis using Qscan 11811 . All 
of the triggers appear consistent with the background popula- 



12 



10 - 




• Foreground 

Background 

RMS Spread 



! tHill 



. . i . i . i i i i i i i i . . . i i i i i i i . i . i . i . i I M 
2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 

CorrPower r 



10 2 



>. 

O) 

i_ 

<D 
C 
HI 

> 
o 

< 

<u 10 
E 



Q. 



c 

> 

LU 



ur 



15 



• Foreground 
— Background 
RMS Spread 




20 25 30 35 

Livingston Normalized Energy Z 



40 



FIG. 9: CorrPower V distribution for background (normalized to the 
livetime of the analysis) and foreground distributions before the final 
CorrPower cut. The gray region is the RMS spread of counts in the 
background time-shifts while the error bars are the error in the mean 
counts per time-shift. The dotted line shows the cut at T=6.2. 



tion. In most cases the triggers arise from the correlation of a 
fairly loud trigger with what appears to be one of a population 
of glitches of smaller magnitude in the other interferometers. 
While only triggers passing category 3 data quality cuts were 
used to set the upper limit, the two events with the highest T 
in the "full" data set after category 2 cuts were also present 
after category 3. Since no triggers in the full data set were in 
apparent excess of the stated upper limits, further follow-ups 
were not necessary. 

In addition to the previously described search requiring data 
from all 3 LIGO interferometers, we also performed a check 
for interesting events during times in which HI and H2 sci- 
ence quality data were available, but LI data was not. The 
two-detector search is less sensitive than the three-detector 
search and background estimation is less reliable, so we do 
not use this data when setting upper limits. However, in 
the first calendar year of S5, there are 77.2 days of livetime 
with only HI and H2 data available (roughly half the livetime 
with simultaneous data from all three interferometers), so it 
is worth checking this data for potential gravitational wave 
candidates. This check used procedures similar to the analy- 
sis previously described, including identical data quality and 
veto procedures. 

Due to the presence of correlated transients in H1H2 data, 
performing time-shifts of one detector relative to the other is 
not a reliable means of obtaining an accurate background. In- 
stead, we use the unshifted H1H2 coincident triggers from the 
H1H2L1 analysis as our estimate of the background since we 
have already determined that there are no gravitational wave 
candidates in this data set. However, the H1H2L1 data set 
is only about twice the livetime of the HlH2-only data set, 
so we are required to extrapolate the false alarm probabil- 
ity distribution to obtained the desired false alarm rate. To 
compensate for the uncertainties in our estimate of the false 
alarm probability introduced by the reduced data set and the 
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FIG. 10: QPipeline significance distribution for background (nor- 
malized to the livetime of the analysis) and foreground distributions 
at the Livingston (top) and Hanford (bottom) sites before the final 
CorrPower cut. The gray region is the RMS spread of counts in the 
background time-shifts while the error bars are the error in the mean 
counts per time-shift. 



extrapolation, we target a more conservative false alarm prob- 
ability of ^0.01 triggers for the HlH2-only analysis. This 
lower false alarm probability and the lack of LI coincidence 
as a veto requires stricter cuts, specifically coherent energy 
Z > 100 from QPipeline and T > 10.1 from CorrPower. 
As in the three-detector search, there were no events above 
threshold (see Fig. fTTT i upon examination of the zero-lag fore- 
ground data, thus no potential gravitational wave candidates 
were identified in the two-detector search. 



VII. SUMMARY AND FUTURE DIRECTIONS 

We have searched the few-kilohertz frequency regime for 
gravitational wave signals using the first calendar year of 
LIGO's fifth science run. No gravitational wave events were 
identified, and we have placed upper limits on the emission of 
gravitational waves in this frequency regime. 

The second calendar year of S5 remains to be analyzed in 
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FIG. 1 1 : CorrPower Gamma zero-lag distribution for the H1H2 anal- 
ysis. The dotted line shows the cut at T=10. 1. 



this frequency range. Several months of this run overlap with 
the first science run of the Virgo II 1 '711 detector, which began 
on May 18th, 2007. During this period of overlap, data from 
Virgo as well as the LIGO interferometers will be incorpo- 
rated into high frequency analysis. Since Virgo is not co- 
aligned with the LIGO detectors, this will require fully coher- 
ent analysis tools rather than CorrPower. Above 1 kHz Virgo 
and LIGO have comparable sensitivities, making their combi- 
nation especially advantageous in the few-kilohertz regime. 

The next LIGO science run will be done with Enhanced 
LIGO Dill , an improved version of the detectors. Most rel- 
evant to high frequency analysis, the dominant background 
of shot noise will be reduced by increasing the power of the 
laser from 10 W to ^35 W, substantially improving the sen- 



sitivity of the detectors. Virgo+, a similarly enhanced ver- 
sion of Virgo, will operate simultaneously. After this, further 
improvements will lead to the AdvancedLIGO l32ll and Ad- 
vancedVirgo 13311 detectors coming online around 2014. Ex- 
tending the analysis of gravitational wave data into the few- 
kilohertz regime will continue to be of scientific interest as 
these detectors become more and more sensitive. 
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